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AllSTUACT 

Observat ions of vibration rotation transitions of silicon monoxide 
in VV f'Ma show that tin* linos originate in accelerating, expanding, 
aiut cool (tUHl 'K) layers of a circumstellar cloud at a distance of 0. 1 l> 
from the central star, the central stellar velocity, as estimated from 
observed SiO 1’ Cypui line profiles, is somewhat redshifted from tin* 
midpoint of the maser emission features. Most of the silicon is prob- 
ably in the form of dust prains. The isotopic ratios of silicon are 
nearly terrestrial. 
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1. IN I'UOPIK no N 

The dense c loud of pas .uul dust surrounding VY Oauis Majoris 
is out’ of the bright ost sources of m.isrr omission aiul infra nxl con 
titmum radiation in our Oalaxy. the complex spatial aiut velocity 
structure of {lie maser emission from Oil, 11. ,0 and SiO has been 
studied by numerous authors (see e. p. ltosen et al. 11)73 and refer 
cnees therein). The optical radiation from \ Y OMa, which is 
scattered by the surrounding nebula, consists of a photospherie, M- 
type spectrum, together with a number of low excitation absorption 
and emission lines which originate outside the photosphere (see e. jj. 
Wallerstein 1071; Herbie, 1 07 4). The optically derived radial velo- 
city of the central star has been uncertain, and has appeared to 
change frequently by larp,e amounts (Wallerstein 1077). This is pre- 
sumably due to the Popple r shiftinp of the photo: pheric spectrum by 
scattering from different parts ot the mo\inp dust shell. Recently 
lteid and Dickinson (197b) deduced a stellar velocity of * 17. 0 l. 9 km s 
(HSR) from the centroid of what is believed to be the thermal component 
of the v 0, .1 3*1 millimeter line of SiO (lhihl et al. 1973) cominp 

from the circumstellar cloud. This value lies approximately midway 
between the OH 1013 Mil/ maser features, and if correct, pives \ Y OMa 
the same velocity as the cluster NOO 3009, to which it may belong 
(llerbip 1999; Humphries 1979), and which is 1.9 Kpc distant. 

In the infrared, low resolution observations between 3.9 i in and 
I t nm (Oillett, Stein, and Solomon 1970) reveal a tl d .0 rather 
featureless spectrum, showinp, that the dust must be present over a 
wide* ranpe of temperatures than in most M piants or superpiants. At 
hiplier resolution near 3.3 in the stellar photosphere is still observed 
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via the overtone bandheads of carbon monoxide (Hyland et al. 1969). 
However, the low excitation and blueshifted velocity of the observed 
fundamental band absorption lines of CO at 4.7 t in (Geballe, Wollman, 
and Rank 1973) indicate that by this wavelength only circumstellar 
material can be seen. 

In order to learn more about the cloud around VY CMa, observa- 
tions have been made near 8.3 Mm of the fundamental vibration rota- 
tion band of silicon monoxide. Although this band is quite complex in 
stellar atmospheres, the cool temperatures and low column densities 
in the circumstellar cloud of VY CMa permit only a relatively few 
lines to be prominent and have allowed a simplified analysis. 


II. OBSERVATIONS 

Five frequency intervals of good atmospheric transmission which 
contain lines of the fundamental band of SiO were observed. The work 
was carried out at the Las Campanas Observatory 2. 5 m du Pont tele- 
scope between 17 January and 26 January 1978. The spectrometer con- 
sists of a liquid nitrogen-cooled Fabry-Perot interferometer with one 
order scanned across the bandpass of a liquid helium -cooled tunable 
grating. The beam diameter of 7", combined with the exit aperture, 
produced a grating bandpass of ~3. 5 cm \ The Fabry-Perot mirror 
spacing was set to give a resolution of 0. 09 cm * (~23 km s *) at SiO 
wavelengths. The lines in VY CMa are partially resolved at this setting. 

Figure 1 shows the spectrum of VY CMa in one of the observed 
intervals, 'the strongest line, which has a P Cygni shape, is from (lie 
ground vibrational state of ^'Si^O. Weaker lines of ^Si^O, 

30 16 p q I e 

Si O and the 2-1 band of ^ Si 1 O are also apparent. Lines of 

28 17 28 18 

"Si O and Si O were not detected and indeed are not expected, 
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In tot.il, sown linos of tho ‘ ' Si O l 0 luiul .uul somowhat fowor 


linos of tin' othor dotootod Minis worn obsorwd bolwoon 1 U’7 oin ' .uul 
123b om ' (S. I inu .uul S. b ‘in). I'ho rot.ition.il lo\ols obsorwd ranp.od 
from .1 2 to .1 33. 


111. ANALYSIS 

Pho moasurod absorption oipiivalont wnltlis of tho SiO linos wrro 
niatohod to isothorm.il slu’ll lino strom.ths to dorivo tho oxoitation 
tomporaturos (for oaoh baud), uolumn ilon.sity, .uul isotopio ratios shown 
m Pablo l. Wlulo it is doubtful that tho linos of oaoh b.unl of SiO ariso 
in layors dosoribablo by a sinp.lo tomporaturo, thorn is insutlioiont 
information on tho donsity ami tomporaluro strnoturo of tho oiroumstol 
lar oloiul to warrant tho uso of a nioro dotailod nnxlol. In addition, tho 
absorption lino shapo is not known ami may bo quito poouliar. 

In tho prosont modol sowral of tho ' SiO l 0 linos aro snlurntod, 


w 


ith optioal dopths noar 2. Pho oolumn donsity m Pablo l was osti 


niatml us in . tho transition probabilitios oaloulatml by Ih'.loUiml ami 

29 . 30 . 

I ambort (1072). IVrivation of tho isotopio ratio Si/ Si is stratp.ht 
forward sinoo linos of both s.pooios of SiO aro of marly m\u.il stronp.th. 

•>0 *' l ) 

Si/ * ' Si was dolorminud by throo mothods, uono ot whioli alono is 

foolproof. In ordor to faoilitato tho oaloulat ions, "oorrootod oipiiva 

lout widths of tho woakost ‘ SiO l 0 absorption linos (whioli aro 

asymnu'trio dm' to noarby omission) wort' oomputod b\ S. . :>linp, tho 

oiiuivalont widths of tho his.h fro.piom y sidos ot thoso linos. Hy ooni 

29 , 

par inn, thoso oipiivalont widths with thoso ot SiO ono obtains 

"^Si, »Si 17 * l. If it is a suimod that own tho hip, h froipionoy sidos 


l 1 
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28 9ft 

of the SiO lines are partially filled by emission, and the /o Si0 absorp- 
lions are widened by giving tliem the same central velocity as '‘SiO 
(see Table 2), one obtains the ratio 23 i 5. Finally, isotopic lines 

no 

and tiie weak 2-1 band lines of ‘ SiO were compared, under the 

p O pQ 

assumption that vibrational LTK holds. Values of ' Si/ Si were 
obtained ranging from 10 to 25 as the temperature varies from 700° K 
to 500°K. A reasonable conclusion from the three results is that 
28 Si/ 29 Si = 20 t 5. 

Velocities of maximum absorption and emission for each SiO 
band are listed in Table 2. They were determined by comparing 
measured frequencies to those calculated from molecular constants 
kindly furnished by D.N.B. Hall (1974). 'Hie frequencies were mea- 
sured relative to observed atmospheric absorption lines and NH^ 
lines. In order to verify that the different absorption velocities found 
for different bands of SiO in VY CMa are real, the molecular constants 
as well as the calibration procedure were checked by calculating the 
radial velocity of a. Orionis from SiO spectra of it obtained during (he 

same observing period as VY CMa. 'Hie l 0 band and 2-1 band lines 

28 +3-1 

of " SiO independently gave the same radial velocity of 23 - km s 

(heliocentric) for a Ori, which is in satisfactory agreement with the 


photospheric value. 

IV. LOCATION OF THE SiO AND 


©F PO° R 


, pAob l- 
OUM-VVY 


PHYSICAL CONDITIONS IN THE CLOUD 
'Hie SiO excitation temperature^ which are near GOO K, imply that 
Uie 8 inn spectrum is formed well beyond the photosphere of VY CMa. 
Distances from the central star of between 0. 13'' and 0. 18" for the SiO 
are obtained from Figure 2 of Ilerbig's (1970) gray body model of (lie 
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oloud and thi’ range of SiO execution (omporaturos m Tablo l of this 
papor. A soeond ilistanoo estimate, whioh usos tin* rotational tempora 
tun* aiul i’stimati'il flux ilonsity at tin* doeonvoKod bottoms of saturalod 
* ‘*SiO linos (using till’ eontinuum flux ilonsity givon by fultolt ot al li> * 0), 
givos a value of 0. lit*, riiis valui* wouKl bo rodueod if a sigmlirnnt 
fraotlou of tho oontimmin radiation noar 8 mu rotnos from tho ilust oat 
shlo of tho SiO. tho inoasuromont of O.lll O.lll* by Mol'arthy ot al 
(1077) for tho radius of VY I’Ma at 8. it in (assuming, a uniform oiroulnr 
energy distr ibution) suggests that this may bo tho oaso. However, a 
radius as largo as O.lll’ would rosult in obsorvod tomporaturos at 
8.H nm noar or bolow 100 K, whioh is s aiowhat lowor than tlio tompora 
tun's found horo. In oonolusiou, most of tho SiO obsorvod noar 8 • in 
probably is looatod noar a distauoo of 0. I ‘a trom tho ooutral star, or 
just outside tho region of radius 0. I whoro H.,H masor emission is 
found (Koseu ot al. 1078). I'ho SiO masors, whioh ro.piiro oxoitation 
lomporatu ros of up to oaOO K (I'uhl ot al. 107 t), probably originate 
olosor to tho star. 

Tho abuudauoo of SiO relative to hydrogen may bo estimated 

15 

roughly from tho oolumn donsity in r.ibU* l, a distance ot l \ U) om 

•)U 

ovor whioh tho ' SiO linos aro formod (assuming a distanoo t v > \ Y t'Ma 

ft 3 

of l. ‘a Upo) and an 11. , density of t x 10 1 om ' . I'ho ilonsity ussuiuos 

u 

0 It 

tho estimated donsity of ill* om in tho region of H.,0 masor omission 
(Kosi'n ot al. 1878) and a l/r*" dopoudonoo. flio rosult is i SiO 11 1 
l \ lit , roughly alt titnos lowor than tho solar valuo tor Si I H !. 
Although this rosult nuiy not bo very accurate, it does suggest that most 
of tho silioon is in dust {’.rains, a condition expected at t h ilistanoo 


from tho ooutral star. 
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The present observations rulo out SiO ns 1 ho source of a broad 
depression between 7. G mil and 9. 3 gm observed in a low -resolution 
spectrum of VY CMa obtained by Gillctt, et al (1970). The feature 
seen by (hose authors has a depth of perhaps 30^, whereas the observed 
SiO band on the average depresses (lie continuum by only In 

addition, the low temperatures found here indicate that (he band is not 
as broad as the depression seen at low resolution. The apparent 
depression may be a result of the emission spectrum characteristic 
of the circumstellar grains. Such emission features have been observed 
in several infrared sources (Russell, Soifer, and Willner 1977; see also 
Allamandola and Norman 1978). 


28, ; 


V. P CYGNI LINES AND THE STRUCTURE AND 
MOTION OF THE CLOUD 

P Cygni line profiles were observed for at least five '“’SiO l 0 

band transistions. The ('mission features must arise in portions of the 

cloud not in (he line of sight to (he continuum source. The emission is 

always redshifted from (he absorption (see Fig. 1), corresponding to an 

expanding cloud. Indeed, the line shape rules out rotation as a large 

component of (he cloud motion in the region of SiO line formation. 

The peak intensities and equivalent widths < i (he emission features 

generally appear much smaller than (he corresponding absorption depths 

and equivalent widths. However, simple models of P Cygni profiles 

observed at the present spectral resolution show that (he deconvolved 

emission and absorption equivalent widths are more nearly equal, with 

the ratio of emission to absorption as much as twice the observed 

28 

ratio. For the few SiO P Cygni lines observed, an inverse relation 


between the ratio of emission to absorption strength and the rotational 



0 


level appears to exist. For example the omission to absorption ratio 

is about one sixth for the l’ 17 lino in Fie,. 1, but about out* half for the 

113 lino (which has a weaker absorption component than Joes the IM7 

lira'). I’he inverse correlation shouUl be confirmed by observing more 

lines. If the inverse correlation were between emission strength aiul 

•)u 

absorption strength, which is also consistent with the ‘ SiO l i) band 
observations, then om' would expect to observe etui- ion features on 


the weak tines of the rare isotopes this was not seen, ('he weak 

% 

' ‘SiO 2 l band also does not show emission, but this may be expected 

since its absorption velocity and excitation indie, tc that it originates 

closer to the exciting source than do the other bands. 

While in principle several mechanisms could lead to the observed 

28 

emission from the first excited vibrational slate lines of ’ ‘ SiO, the 
dominant one in this case is simply reemission (resonant scalterine,). 

Hy usine, equations {*,iven by Millikan and While (1983) for computing; 
collisonal relaxation times, it can be shown that at the densities expec 

i) 

toil in the emitting part of the cloud ( U)‘ cm ) the collisional vibra 
tional excitation rate is only 10 ' of the radiative excitation rate at 
8 am. Vibrational excitation by \ am or shorter wavelength photons 
followed by relaxation \ ia 8 -m emission is also relatively unlikely be 
cause of the low transition probabilities for overtone band absorption. 

For resonant scattering in a spherically symmetric and expanding 
cloud without dust, symmetry and conserv 'alien of energy require that 
the line absorption and omission observed from any lireetion be equal, 
if the central continuum source is small. If, as may well be the ease 
for VY t’Ma at 8 am, the continuum source is not small compared to 
the SiO cloud, the emission is reduced duo to the absorption of scat 
ton'd radiation by the continuum source. Oust in the line forming 
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region absorbs photons reemittcd by SiO, and further reduces the 

emission, most strongly for the optically deepest transitions. The 

limited data for the SiO 1-0 band is in agreement with this model. 

However, the lack of emission for the optically thin transitions of 
29 „ 30 

SiO and SiO is puzzling. It is possible that much of the absorp- 
tion equivalent widths of these lines may originate in layers which are 

close enough to the opaque surface that the emission is highly sup- 

28 

pressed. In contrast SiO may absorb a substantial fractional 
amount of radiation in cooler outer layers, where reemitted photons 
are more likely to escape. 

Clearly, a detailed model and observations at higher spectral 
resolution are necessary to better understand the observed line 
shapes. However, on the basis of the present observations, it is 
possible to draw some conclusions about the shape of (lie circumstellar 
cloud and its central velocity. The presence of some emission fea- 
tures at least one-half as strong as the corresponding absorption 
lines limits the possible geometries of the expanding cloud. Some 
authors (Herbig 1969; van Blerkom and Auer 1976) have proposed 
that the circumstellar material at certain distances from the central 


star is in the form of a disk seen nearly edge on. The resonant 
scattering mechanism, would, in the case of a disk of this orientation, 
produce in our direction a very weak SiO emission line, whose maxi- 
mum (small continuum source, dust free) equivalent width relative to 
the associated absorption would be roughly the ratio of the thickness 


of the disk to its diameter. Thus, the maximum observed ratio of 

28 

— 1 / 2 rules out a thin, nearly edge-on disk, near the radius of * SiO 


line formation (r ~ 0. 15'). 
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Hie central stellar velocity may 1 c iMimalcd from the 1’ Cygni 
profiles in the following way. For the case of a cloud expanding sym- 
metrically about a continuum source and with the exceptions noted 
below, the central velocity of the emission feature corresponds to the 
velocity of that source. Fits to the observed profiles based on various 
simple models of P Cygni line formation show that the observed velo- 
city of peak emission is redshifted by 7 ' 3 km s * from the central 
velocity of the emission feature (due to the present instrumental reso- 
lution). Thus, the observed velocity of peak emission (39 » 4 km s *) 
implies a central velocity of 32 1 7 km s * (LSP). 

Two effects not included in the simple models could alter this 

result. Iteabsorption by SiO in outer parts of an accelerating cloud (see 

% 

below) will decrease the apparent emission from the front half causing 
the peak emission to occur behind the limb and be mure redshifted 
than the central velocity. Absorbing dust within the line forming region 
will have the opposite effect, reducing the contribution from the rear 
parts of the shell more than from the front. Scattering of radiation by 
dust external to the SiO could shift the 8 am spectrum to the red in a 
manner similar to that proposed by llerbig (1989) and calculated by 
van Blerkom and van Hlerkom (1978) for the optical spectrum of VY I’Ma. 
However, the scattering cross section of normal si/.e dust grains is 
negligible at 8am. To the extent that the first two effects are small or 
compensate, and for the cloud shapes which are either spheres or 
disks, the above velocity should represent (ho slell ir velocity. 

Two lower limits to the stellar velocity can be obtained from the 
8 am data and are consistent with the al we value. The midpoint of the 
entire P Cygni profile (absorption and emission), averaged over seve- 
ral lines, is 25 t5 km s . This value is a lower limit to the central 
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velocity because emission from the most redshifted SiO is blocked by 

the continuum source. The second lower limit is (he most redshifted 

“1 

absorption velocity, 18 * A km s which is that of the ' SiO 2-1 band. 

Since the SiO lines are formed far from th- stellar photosphere, the 
gas containing tlie observed SiO is expected to be moving away from 
the star. This implies that the stellar velocity is redshifted with 
respect to the velocity of the 2-1 band. 

The above arguments indicate that the stellar velocity, although 
in between the 1612 MIIz OH features, is probably not central to them. 
Only the lower of the two lower limits, from the 2-1 SiO band, is 
consistent with l?eid and Dickinson’s (1976) proposed stellar velocity 
of 17. 6 il. 5 km s \ which is central with respect to the maser 
velocity structures. The stellar velocity of ~32 km s proposed here 
apparently requires some asymmetry in the expanding envelope, such 
that no 1612 MHz OH emission from the extreme rear of the envelope 
is directed towards us at present. In addition, it requires that most 
of the H^O and SiO maser emission arises in the front half of the cloud. 
It is interesting that the absorption line velocities of the various 

bands of SiO are not all the same (see Table 2). The velocity of the 

28 - 1 
ground vibrational state lines of SiO is blueshifted by about 13 km s 

relative to that of the first vibrational slate lines. 'Hie latter lines 

must arise in a region interior to the former lines; thus, the velocity 

difference corresponds to outward acceleration. This interpretation 

is consistent with that of the P Cygni lines. The intermediate velocity 

of the SiO and ’ SiO lines indicates that they form largely in between 

the two "°SiO bands. The expansion velocity at r ~0. i 5* where 8 urn 

SiO lines originate is roughly half the expansion velocity at r ~l" where 

OH 1612 MHz emission originates. 
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VI. THE SILICON ISOTOPIC HATIOS 
28 29 28 30 

The terrestrial ratios ‘ Si/ ' ' Si and Si/ Si are 20 and 1.5, 
respectively. Within the uncertainties, the results for \ V CMa differ 
from these only in a probable slight enhancement of '^Si. Silicon 
isotopic ratios are available at present in only two other objects out- 
side of the solar system (Peer, Lambert, and Sncden 1974; Clark and 
Lovas 1077) and are also nearly terrestrial. Hie re nits an begin- 
ning to suggest a uniformity of these ratios over a large region of our 
Galaxy. 

Isotopic abundances of silicon in a stellar interior are not expected 
to change until advanced stages of red giant stellar evolution are 
reached. Only as early as helium shell burning might s process 

reactions change (increase) the abundances of the rare isoto; \s rcla 

28 

tive to " Si. In the still later stages of oxygen burning and silicon 

28 29 30 

burning, cores are produced in which ‘ ' Si/ ’ Si and Si/ Si are 

larger than their terrestrial values (Podansky, Cla>ton aiul i'owler 

\ 

1968, Woosley, Arnett and Clayton 1972). All of these mechanisms 
require some mixing process to alter the surface abuiulane.es. 

Whether VY CMa is a newly formed star or an evolved object has 

% 

been a matter of considerable discussion (e. g. see Herbig 1969; Hyland 
et al. 1970; Geballo et al. 1978). The observed isotopic abundances of 
silicon might restrict VY CMa to an evolutionary phase earlier than He 
shell burning, but certainly do not rule out earlier red giant phases. 
'Puis, the present results in themselves leave open thee stion of the 
evolutionary state of VY CMa. 
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SUMMARY 

1. The 8 M* in SiO lines in VY CMn arc formed in ~G00°K accelerating 
and expanding gas at a radius of -0. 15 ". In this region, the gas 
is not in the shape of a thin disk seen nearly edge -on. 

2. The radial velocity of the central star appears to be shifted by 
about r 15 km s * from the midpoint of the 1612 Mil/ OH pattern. 

3. Most of the silicon probably is in the form of dust grains. 

4. The silicon isotopic ratios are roughly terrestrial. 
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